Evaporatively driven morphological instability 



(N 

o 

(N 

a 

(N 



o 

£ 



I 

c 

o 
o 



> 

(N 

m 
cn 

o 



X 



Robert W. Style 1 and J. S. Wettlaufer 2 
1 Institute of Theoretical Geophysics, Department of Applied Mathematics and Theoretical Physics, 
University of Cambridge, Wilberforce Road, Cambridge, CB3 OWA, U.K. and 
2 Department of Geology & Geophysics and Department of Physics, 
Yale University, New Haven, Connecticut 06520-8109, USA 
(Dated: 23 June, 2007) 

Simple observations of evaporating solutions reveal a complex hierarchy of spatio-temporal insta- 
bilities. We analyze one such instability suggested by the qualitative observations of Du and Stone 
and find that it is driven by a novel variant of the classical morphological instability in alloy solid- 
ification. In the latter case a moving solid-liquid interface is accompanied by a solutally enriched 
boundary layer that is thermodynamically metastable due to constitutional supercooling. Here, we 
consider the evaporation of an impure film adjacent to a solid composed of the nonvolatile species. In 
this case, constitutional supercooling within the film is created by evaporation at the solution-vapor 
interface and this drives the corrugation of the solid-solution interface across the thickness of the 
film. The principal points of this simple theoretical study are to suggest an instability mechanism 
that is likely operative across a broad range of technological and natural systems and to focus future 
quantitative experimental searches. 



I. INTRODUCTION AND MOTIVATION 

Evaporation of water is an important phase transfor- 
mation appearing in many guises throughout everyday 
life. Its effects range from processes intrinsically linked 
to the Earth's hydrologic cycle, to the regulation of body 
temperature in hot environments, to the production of 
coffee-ring stains beneath a spilt coffee droplet and to 
important processes underlying microfluidics. While the 
evaporation of bulk pure materials provides a wide vari- 
ety of thin film instabilities [l| , the simple act of introduc- 
ing solutes, with their myriad of interfacial effects, leads 
to a spate of new behaviors whose exploration is in its in- 
fancy. Examples include the investigation of coffee-ring 
formation by a sessile, particle-laden droplet Q , the ob- 
servation of finite contact angles in evaporating, wetting 
films [3, 3 , and of particular interest here, the observa- 
tions of Du & Stone on evaporatively grown salt trees Q . 
Understanding the underlying mechanisms that control 
the patterns that emerge during evaporation and solidifi- 
cation of solutions influences a wide range of fields, from 
the interpretation of evaporite deposits for clues about 
environmental conditions in the geologic past [f|, to the 
scaling deposits in tea kettles when used with hard water 
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Du & Stone [f| observed the evolution of ramified 
crystal structures during the evaporation of water from 
an ammonium chloride solution under varying (but not 
quantitatively controlled) external conditions. Initially 
these crystals were confined within the bulk liquid, but 
eventually they protruded above the liquid surface, form- 
ing salt trees, the growth in this latter stage being caused 
by evaporation of a film maintained at the surface of 
the solid salt by capillarity. In consequence Du & Stone 
@ postulated the existence of two distinct instabilities 
giving rise to two different characteristic lengthscales of 
the system (figure H]). Although their observations did 
not quantify the conditions for onset of either instabil- 



ity they were led to make the following general distinc- 
tions. The smaller scale instability (with a characteristic 
length scale estimated to be ss 10/im) was suggested to 
arise from evaporatively driven supersaturation that al- 
lows fluctuations of the salt-liquid interface to grow in 
a manner qualitatively similar to the growth of a solid 
into a supercooled melt. The second instability (with a 
characteristic length scale estimated to be ss 100/^m) was 
proposed as being due to a spatial variation in evapora- 
tive flux across the surface; the film at the surface of a 
crystal protruding further into the dry atmosphere will 
cause faster evaporation of the film and hence further 
growth relative to adjacent regions. 

While these rough estimates neither constrain the na- 
ture of either instability, nor quantify the thermodynamic 
ranges of their existence, the Du & Stone @ observations 
motivate a simple theoretical analysis that we hope will 
focus the future experimental search. Indeed, more so- 
phisticated treatments are not warranted by the present 
experimental evidence. Hence, in this paper we focus on 
the first of these instabilities and demonstrate the exis- 
tence of an observable length-scale during the evaporative 
growth of salt crystals. 



II. THEORY 

We study the mechanism of the first, salt-diffusion- 
driven, instability and simplify the system as follows. 
Firstly, although latent heat is produced at two phase 
boundaries, because the diffusivity of heat in the liquid 
phase much greater than the diffusivity of salt D, 

we arc justified in treating the system as isothermal. Sec- 
ondly, the vapor pressure of salt in air is exceedingly small 
relative to that of water vapor and hence there is no salt 
transport in the air. Additionally, apart from the conse- 
quences of the instability we discuss, there are no liquid 
inclusions in the salt phase. Thirdly, we treat the system 
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FIG. 1: Scanning electron micrographs of typical ammonium 
chloride salt trees, (a) A branch of a dried tree approxi- 
mately 3mm in length. The picture demonstrates the multi- 
scale roughening resulting from the evaporative process, (b) 
Cross section of the porous interior of a typical dried tree. 
Roughening is caused by evaporation in the presence of a 
wetting film [f| . These images motivate our study but do not 
provide a length scale against which theory can be quantita- 
tively compared. 



as planar with shallow perturbations to the two inter- 
faces. This allows us to use the typical approximations 
to the curvature and normal derivatives that we check 
post hoc. Finally, we assume that the liquid-vapor inter- 
face is maintained flat by flow in the liquid phase due to 
the influence of surface tension. From lubrication theory, 
we find that the time scale for the relaxation of a film of 
thickness do to a planar state is given by 



(1) 



where p is the dynamic viscosity of water, a is the 
wavenumber of the perturbation, and o~i a is the surface 
tension of the liquid-air interface. The diffusive time 
scale for the system is given by td = d\jD where D 
is the diffusivity of salt in the liquid. Thus the hydro- 
dynamic relaxation to a planar interface will occur on 
a much shorter timescale than that on which instability 
occurs as long as 



< 1, 



(2) 



where uj max is the maximum growth rate of the insta- 
bility as given later. Hence, because a ~ 27r/10 _5 m (as 
we show later) , we can ignore the motion of the liquid- 
vapor interface when do 3> 6 x I0~ lo m. This is a good 
assumption in the relevant parameter space. We will also 
show that the advection of solute in the film is negligible 
relative to diffusive transport, and so the precise details 
of the flow can be ignored. These are either common 
approximations that are borne out experimentally or are 
confirmed to hold within the regimes of validity of the 
theory. Finally, although we have analysed the effect of 
marangoni stresses in transporting liquid to the edges 
of evaporating droplets @ previously, the salt concen- 
tration at the liquid-vapor interface was taken as fixed, 
giving rise to large salinity differences across the liquid 
film. Here, we relax this assumption and consider instead 
the evaporatively controlled salinity at the liquid-vapor 
interface. Motivated by the observations of Du & Stone 
we shall treat the ammonium chloride-water system 
at room temperature and typical humidities throughout 
this paper. However, we note that the analysis is appli- 
cable to many other systems of the same general class. 

Consider a planar salt crystal covered by a film of salt 
solution that has an initial thickness do as shown in fig- 
ure @. We expose the film to an undersaturated vapor 
which drives evaporation at the liquid-vapor interface. In 
general, within the film the solute density then satisfies 
the advection-diffusion equation 



dp ~ 



D\7 2 p 



(3) 



where p(x, z, t) is the density of solute in the liquid, u is 
the velocity of the liquid and x, z, t are horizontal, vertical 
and time coordinates. 

Two-phase, two-component equilibrium at the two 
phase boundaries, requires the satisfaction of two bound- 
ary conditions for the system. Firstly, the liquidus rela- 
tionship, which defines bulk equilibrium between the salt 
in the film and that in the underlying planar crystal, is 
modified for the shift in equilibrium associated with the 
curvature of the solid-liquid interface 1C s i as 



T = T + mp - 



C 



-K 



si •■ 



(4) 



where T m is the melting temperature of a pure liquid salt, 
m is the slope of the liquidus, a s i is the surface energy 
of the salt-liquid interface and C / is the heat of solution 
of the salt. This describes equilibrium as a linear combi- 
nation of the colligative and Gibbs- Thomson effects. We 
note that we have assumed that the curvature term is 
independent of solute density, so the resulting expression 
will not be exact. Rearranging, we obtain 



T — T 



m 



0~ sl^rn p. 

mCf 



(•5) 



which gives us the equilibrium boundary condition at the 
salt-liquid interface: 

p{h s i) = p° B + jiC s i. (6) 
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FIG. 2: Schematic for the planar growth of a salt crystal due 
to the evaporation of a salty film, p~B and pr are the den- 
sities of salt in the solution at the salt-liquid interface and 
the liquid-vapour interface respectively, E is the evaporation 
rate of the film, and Uo is the growth rate of the salt crystal. 
Initially the interfacial positions are given by z = for the 
salt-liquid interface and z — do for the liquid-vapour inter- 
face. 



The natural length, velocity and time scales imposed 
on the system are given by do, E and do/E respec- 
tively and so we scale variables accordingly so that 
(x, z, hiv, h s i) = d fa, z, hv,h s i), (u, U ) = E(u, U ) and 
i = dot/E. Thus, equation ([3]) becomes 



u.Vp ) = V 2 p, 



(10) 



where the Peclet number Pe is defined as d$E/D and 
corresponds to the relative importance of advective (due 
to evaporation of the interface) and diffusive transport 
processes. For a typical aqueous ammonium chloride so- 
lution at 20°C, E ~ lO^ms" 1 and D = lO^Vs" 1 
so that Pe~ do(m) x 10 2 - Therefore, the Peclet number 
will be small for all films under consideration and we can 
approximate equation (J3J by Laplace's equation 



V 2 p 



0. 



(11) 



Hence, during planar growth of the salt crystal informa- 
tion is diffused sufficiently rapidly across the film that 
the solutal profile is linear and given by 



(P°T 



(12) 



Here, p^ is the solute density at the liquid- vapor inter- 
face in this basic state, determined from the boundary 
conditions ([SJ and ^ upon application of which we find 



h s i is the vertical position of the solid-liquid interface 
relative to the position of the planar growth interface 
position as derived below, p B = (T m — T)/m is the equi- 
librium value of the solute density at a planar interface, 
and 7 = a s iT m /mCf « 1.8 x 10~ 7 kgm~ 2 . 

Secondly, the dynamically maintained planarity of the 
liquid-vapor interface allows us to treat the evaporation 
rate there as a constant. Although the evaporation rate 
may be influenced by gradients in salinity along the sur- 
face, these are shown below to be small because p « p B . 
Thus we take 



hi v — E. 



(7) 



The two boundary conditions needed to close the sys- 
tem of equations express conservation of salt across the 
solid-liquid interface, 



D 



dp 

dz 



{h lv + U )p s 

all j 



(8) 



and conservation of water across the liquid-vapor inter- 
face, 



az 



E P\z=h, 



(9) 



Here p sa it is the density of salt in the solid phase and 
Uq is the planar growth rate of the solid-liquid interface 
before the introduction of any perturbation. 



P°B 



1 



and 



Uop S al 



1 - Pc 



Ep° T . 



(13) 



(14) 



We note that this last equation agrees with our intuition 
expressing as it does the fact that as water is evaporated 
at the surface, the amount of salt that was dissolved in 
the evaporated water is instantaneously deposited onto 
the salt crystal at the base of the film. Also, from equa- 
tion[13]the gradient in salt density across the film is given 
by 



P°T 



Pb 



p^Pe _ p B E 
dn ~ D 



(15) 



The last approximation also stems from the fact that 
the Peclet number is small. This shows that the salt 
density gradient across the film is roughly independent 
of film thickness, and thus any perturbation at the salt- 
liquid interface will encounter similar salt field conditions 
independent of film thickness. Thus, in this limit, it can 
be anticipated that the growth rate of the perturbation 
will not depend on do- 

We are now in a position to conduct a linear stability 
analysis of this quasi-stationary planar state, noting first 
that for this quasi-stationary linear stability analysis to 
be relevant, the perturbation must grow much faster than 
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the rate of change of the underlying base state, so we 
must check a posteriori that the growth rate of the most 
unstable wavelength satisfies a max > p^ I Psalt- From 



consideration of the basic state, we choose to nondimen- 
sionalise the solute density by setting p—p% — (p T ~p B )p. 
The governing equations (|3|6|8|9j) then become 



V 2 P = 0, 



p = VK S 



dp 

dz 



z—h s 



Psalt 



and 



Here, L 



dp 

dz 

7/($ 



Pep| 



z=hi 



Pi 
P°T 



(16) 



(17) 



(18) 



(19) 



j t — p B )do is the nondimensional solutal 
Gibbs-Thomson coefficient, and we recall that the in- 
terface position hi v = 1 — t. We now impose a small 
perturbation upon the interface so that 



h s i = ee lax+at . 



(20) 



The perturbation to the solutal field that satisfies 
Laplace's equation is given by 



p = z + (Asinhaz + £?coshaz)e 



iax+at 



(21) 



and after applying boundary conditions (|17H19|) . we ob- 
tain 



aA 



Psalt 

Pt 



B = e(Ta 2 - 1) 



(22) 



(23) 



and 



A sinh i 



Aa cosh a Ba sinh a 
Bcosha = h - (24) 



Pe 



Pe 



from which we obtain the dispersion relation governing 
the instability: 

Pt ^ n 2\ /atanha - Pe\ 

a = J - J —a{\ - Ta z ) — . (25) 

; 1 - - Petanha ' 



Psalt 



This relation provides the most unstable wavenumber or 
cutoff wavelength as 



1 

7f' 



(26) 



which for a film thickness of do = 10~ 4 m takes a value 
of a co — 45.3. Due to the fact that this is much greater 
than 1, we suggest that a reasonable approximation to 
the dispersion relation is to take the limit for large a, 



TABLE I: Table of typical values for the ammonium chlo- 
ride/water system 



Constant 


Value 


Units 


D 


77^ — TT" 

10 J 


t> r 

nrs 


E 


10 7 


i 

ms 


7 


1.8 x 10~ 7 


1 2 

kgm 


psalt 


1.5 x 10 3 


kgm -3 


P°B 


3.7 x 10 2 


kg m -3 


£f 


4 x 10 s 


Jm" 3 




611 


K 


m 


0.25 


Km 3 kg" 1 




3 x 10~ 2 


Jm" 2 


Ola 


7.6 x 10" 2 


Jm" 2 


P 


1.787 x 10~ 3 


kg m s _1 



in which case it reduces to the same form as that of the 
Mullins-Sekerka instability (e.g., @), namely, 



Psalt 



a(l - La 2 ). 



(27) 



Figure (jTT]) shows a plot of the dispersion relation for 
c?o = l()- 4 m and also the asymptotic approximation 
given by the equation above. It can be seen that the 
asymptotic approximation gives excellent agreement with 
the full expression. Indeed for film thicknesses greater 
than a micron, the asymptotic expression can be substi- 
tuted for the full expression without appreciable loss of 
accuracy. 



<J 




LIG. 3: The dispersion relation for the instability with do = 
10 _4 m. Insert shows the behaviour close to the origin. Con- 
tinuous line gives full dispersion, dashed line shows large a 
asymptotic relation. 



From equation (|27l) , the most unstable wavenumber is 
given by 



(28) 



3L 



and thus we find the growth rate of the instability, 
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a{a r 



to be 



Psalt 



(29) 



Thus for a typical film of thickness do = 10 _4 m, the 
wavelength of the instability will be 2.5 x 10~ 5 m with a 
growth rate of around 16 times the speed of growth of 
the crystal. 

The most unstable wavenumber corresponds to a wave- 
length 



A cx doVT — y/l^do, 



(30) 



where Z 7 = j/(p^ — p%) is the capillary length corre- 
sponding the supersaturation across the film. This ex- 
presses the rule of thumb that the length scale char- 
acteristic of morphological instabilities is the geometric 
mean of the diffusion length and the capillary length (e.g, 

BUS)- it 

is also informative to note the dependence of 
and a m ax upon do- 
In the small Peclet number limit, equation (|13[) be- 
comes 

p° T - p% « Pep^ (31) 
so that the most unstable wavenumber is (dimensionally) 



3jD 



(32) 



and the dimensional growth rate of the instability is 



2E p%E 
3 V 3 7^ ' 



(33) 



Therefore, as was suggested by the fact that the salt gra- 
dient in the film is independent of film thickness, both 
of these results are independent of do and their sole de- 
pendence on the diffusive driving force, the evaporation 
rate E, is simply exhibited. Although Du and Stone [j| 
did not measure the evaporation rate, here we estimate 
a value of 10 -7 ms -1 , from the typical room tempera- 
ture situation in which they made their observations, and 
hence we find a wavelength of 23pm. 



III. CONCLUSION 

We have illustrated the existence of a new instability in 
the evaporation of impure (salty) films. Supersaturation 



of the film caused by evaporation at the liquid-vapor in- 
terface leads to enhanced growth of perturbations to the 
salt-liquid interface underlying the film. The instability 
is similar to instabilities arising during the solidification 
of a binary alloy into a melt [9j , however the chief differ- 
ence is that supersaturation is created at some distance 
from the advancing salt front, whereas with the classical 
morphological instability of directionally solidified alloys, 
the supersaturation is caused by rejection of solvent im- 
mediately adjacent to the salt front. As suggested by 
Du and Stone [j| based on qualitative observations, this 
is a likely mechanism for the growth of the microscopic 
perturbations to the surface of growing ammonium chlo- 
ride trees. The analysis is appropriate for small Peclet 
number which, given typical expected room temperature 
evaporation rates, requires films less than O(lmm) in ini- 
tial thickness. It should also be noted that solutal con- 
vection would be expected to set in for films on the order 
of 1mm in thickness. However we expect films present in 
ammonium chloride trees to be sufficiently thin to avoid 
these complications. 

Although Du and Stone Q observed structures on the 
order of 10//m, our analysis is principally intended to 
point out the nature of this evaporative instability and 
the experimental parameters of note, namely, 7 and E. 
The wide variability in the observed scales across the sur- 
face of the film is likely to be due to gradients in evapora- 
tion rate at the interface due to local conditions and/or 
an underlying nonlinear bifurcation (common in diffusive 
type solidification problems) but the present state of the 
experimental evidence does not warrant a more thorough 
theoretical treatment. Thus, while this mechanism pro- 
vides the likely cause of the growth of microscopic per- 
turbations as observed by Du and Stone 0] , it clearly co- 
exists with the other phenomena that they have pointed 
out and we have discussed further here. It is hoped that 
by focusing on a well defined situation, this analysis will 
provide a framework for refined experimental searches. 
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